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ESTIMATION  OF  RADAR  RANGE  BIASES  AT  MULTIPLE  SITES 


INTRODUCTION 

An  interesting  and  important  problem  the  Navy  has  been  concerned  with  for  years  is  called 
“Gridlock.”  Primarily,  it  is  concerned  with  using  radars  to  form  a  relative  navigation  system.  This 
can  be  achieved  because  radars  located  on  different  ships  should  see  common  targets  at  the  same 
points  in  space.  In  general,  there  is  a  large  number  of  identifiable  parameters  which  must  be  esti¬ 
mated  such  as  latitudes,  longitudes,  speeds,  heights,  orientations,  range  biases,  azimuth  biases,  etc. 
Since  these  parameters  are  often  coupled,  often  the  general  estimation  problem  can  become  very 
complex.  It  would  be  desirable  to  formulate  the  problem  so  that  some  of  the  parameters  could  be 
estimated  independently  from  the  rest.  If  this  could  be  done,  several  small  estimation  problems 
could  be  solved,  and  probably  a  more  stable  solution  could  be  obtained. 

A  means  of  estimating  the  range  bias  error  which  is  independent  of  a  ship  or  platform’s  posi¬ 
tion  and  orientation  is  studied.  The  method  to  be  described  is  based  on  the  requirement  that  two 
non-colocated  radars  must  measure  the  distance  between  two  targets  as  the  same.  The  algorithm  is 
first  formulated.  The  performance  of  the  estimation  is  studied  using  synthetically  generated  data 
and  real  recorded  data.  The  recorded  data  were  obtained  simultaneously  using  SPS-39  radars  located 
at  the  Johns  Hopkins  University/ Applied  Physics  Laboratory  (JHU/APL)  and  the  Naval  Research 
Laboratory  (NRL),  Chesapeake  Bay  Detachment  (CBD).  Finally,  the  results  are  summarized. 


MATHEMATICAL  FORMULATION 

Figure  1  illustrates  the  geometry  of  the  problem.  Coordinate  systems  are  established  at  each 
radar  site  with  the  x-axis  oriented  due  east,  the  y-axis  oriented  due  north  and  the  z-axis  oriented 
outward  along  a  ray  joining  the  site  and  the  Earth’s  center.  Azimuth  is  measured  counterclockwise 
from  the  x-axis,  elevation  is  measured  from  the  x-y  plane  and  range  is  the  length  of  the  position 
vector  whose  initial  point  is  the  origin  of  a  radar  site’s  coordinate  system  to  a  point  in  space. 

The  difference  vectors  from  sites  1  and  2,  Dx  and  D2,  can  be  described  as  functions  of  the 
following  elements: 

•D1  -  j,  rj  2,  1  j ,  aZj  2,  j ,  e®x  2^  * 

and 

D 2  F( r2i < r22’  aZ2\ ’ °*22’ e®21’  e^22^ ’ 


(1) 

(2) 


Manuscript  submitted  May  22,  1981. 


i 


1 


NRL  REPORT  8510 


and 


rT  =  the  true  range  measurement  for  target;';  from  radar  i, 

1  ii 

b  =  the  range  bias  in  radar  i, 

ri 

n  =  the  noise  present  in  radar  i’s  range  measurement  to  target;', 
ru 

azT  =  the  true  azimuth  measurement  for  target ;';  from  radar  i, 

1  a 

b  =  the  azimuth  bias  in  radar  i, 

az  i 

n  =  the  noise  present  in  radar  i’s  azimuth  measurement  to  target ;, 

aZjj 

=  the  true  elevation  measurement  for  target;;  from  radar  i, 

1  a 

bon  =  the  elevation  bias  in  radar  i, 

exi 

and 

n  n  =  the  noise  present  in  radar  i’s  elevation  measurement  to  target;. 
ex0 

Referring  again  to  Fig.  1,  the  location  of  target ;  from  radar  i,  T is  given  in  vector  form  using  the 
unit  vectors  i  , ;  ,  and  k  by 

Ti  i  ~  rlt  cos  eCj  t  cos  aZj  j  i  +  r}  x  cos  eCj  x  sin  az1  j  j  +  rx  1  sin  eCx  x  k  ,  (4) 

and 

^12=  rl2cose^l2cosazl2^*  +  ri2  c°s  e%i2  sin  a*l2T  +  ri2sin  e%!2^  (5) 

as  seen  from  radar  1  and 

?21  =  r21  cos  eS21  cos  az21  7*’  +  r21  cos  e£21  sin  az2X  T  +  r 21  sin  eC21  h  ,  (6) 

and 

T22  =  r22  cos  ei22  cos  az22  i  +  r22  cos  e8,,„  sin  az22  j  +  r22  sin  e£22  k  (7) 

as  seen  from  radar  2.  The  distance  between  targets  1  and  2  as  measured  by  radar  1  is 

D\  =  |T12  "  Tll|2  ’  <8) 
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and  as  measured  by  radar  2  is 


^2  =  1^22  ~  -^21  (  ’  (9) 

where  |  |  represents  the  magnitude  of  the  vector. 

Using  Eqs.  (4)  through  (7),  the  distances  Dj  and  D2  are  given  by 


D1  =  rll  +  r12  ~  2rllr12[sin  e^21  sin  e®22  +  cos  e®21  COS  e®22  cos  (az21  ”  “^22^]  ’  (10) 

and 

D\  =  r21  +  r22  ~  2r21r22  fsin  e821  sin  eC22  +  cos  eg21  COS  eg22  COS  (az21  '  ^22^1  *  (J1) 


Approximations  for  Dx  and  D2  are  obtained  by  expanding  the  functions  in  a  truncated  Taylors 
series  in  terms  of  the  noise  and  bias  errors.  If  the  function  D1  and  D2  are  nearly  linear  functions  in 
terms  of  range,  azimuth,  and  elevation  over  the  range  of  values  of  the  noise  and  bias  errors,  only  the 
first  two  terms  of  the  Taylor  series  need  be  kept  and  the  approximation  becomes  linear.  In  the  case 
we  will  consider,  the  elevation  and  azimuth  bias  errors  are  assumed  to  be  zero  or  previously  removed 
by  another  process.  Actually  it  can  be  shown  that  D1  and  D2  are  not  dependent  on  a  constant 
azimuth  bias.  Assuming  the  elevation  bias  is  zero,  the  linear  approximation  for  D 1  is 


D 


1 


(12) 


where 


- 1 

§ 
t— * 

_ 1 

ry.  T  _ 

'  ailj 

dDl 

O’ 

b 

•—* 

— i 

Lav 

•LC'*J 

dn 

12 

dn 

i  2 

an^12j 

4 
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Dl  is  the  approximate  value,  D}  is  the  true  or  mean  value  at  the  same  time,  is  a  1-by-l  matrix 
derived  in  the  appendix,  and  2  is  a  matrix  of  partial  derivatives  which  is  also  derived  in  the 
appendix.  Similarly,  D2  can  be  approximated  as 


(13) 


where 


r  T 

3D2 

[c„  1  = 

L  21J 

a  br 

r2  m 

3D 2  3 D2  3D2 


3nr  3 n,  3 n 
r21  r22  a 


‘21 


3D2 


If  the  measurements,  Eqs.  (12)  and  (13),  are  made  in  time  coincidence,  the  residue  error 
A D  =  Dx  -  D2  between  the  measurements  is 


Equation  (14)  can  be  rewritten  as 

AD  =  HX  +  N , 


(14) 


(15) 
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where 


and 


"■[[c.«]'[c«]]'  x' 


W'[[C>2]-[C22]] 
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'  22 


>«22 


Kquation  (15),  without  the  noise,  is  the  form  of  one  linear  equation  with  two  unknowns. 

If  another  equation  using  at  least  one  other  target  is  formed  like  (15),  we  would  have  two  linear 
equations  and  two  unknowns  which  would  have  a  unique  solution.  This  system  of  equations  is 
represented  by 


R  =  GX  +  S , 


(16) 


where 


'ad,' 

V 

■/l 

r! 

'< 

R  = 

,G  = 

,x  = 

,  and  S  = 

hk 

>2. 

.N\ 

The  subscripts  on  D,  H,  and  N  denote  the  values  in  Eq.  (15)  corresponding  to  the  first  equation  and 
second  equation  formed  by  examining  distances  between  different  target  pairs.  Equation  (16)  is 
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now  in  ihe  form  of  the  observation  equation  used  in  the  least-square  estimator.  The  vector  of  range 
bias  errors,  .Y,  is  the  state  vector  with  S  being  a  matrix  of  measurement  noises. 

The  smoothed  range  bias  errors  are  estimated  recursively  by  a  least-square  filter  from  Ref.  1 . 
The  three  steps  involved  in  the  recursive  algorithm  are  as  follows: 

Step  1.  Calculate  the  prediction  observation  A (k)  ; 

A(fc  +  1)  =  A(/?)  +  G'C~  lR(k)  ;  (17) 

Step  2.  Calculate  a  new  covariance  matrix  Pi  k  +  1 ), 

P~1(k  +  1)  =  P~l(k)  +  G'C^G  ■ 

Step  3.  Calculate  a  new  smoothed  estimate  x  (k)\ 

X(k)  =  P(k  +  1  )A(/r  +  1) 
where  C  is  the  covariance  matrix  of  the  noise  S. 

The  following  sections  describe  the  performance  of  this  algorithm  with  simulated  and  real  data. 


SIMULATION  RESULTS 

The  simulation  used  to  test  the  algorithm  consisted  of  two  stationary  platforms;  platform  one 
had  (x,  y,  z)  coordinates,  (10,  8,  0)  measured  in  nautical  miles,  while  platform  two  was  located  at 
the  origin.  On  each  simulated  iteration,  two  stationary  targets  were  randomly  generated  in  any  one 
of  the  four  quadrants.  The  radar  at  the  two  platforms  then  made  their  own  range,  azimuth,  and 
elevation  measurements  on  each  target  as  illustrated  in  Fig.  1.  All  of  these  measurements  were 
corrupted  with  zero-mean  noise  in  range,  azimuth,  and  elevation  with  standard  deviations  of  360  m 
(1200  ft),  0.5°,  and  1°  respectively. 

The  filter’s  ability  to  accurately  approximate  range  bias  was  tested  using  the  bias  conditions: 

Platform  1  Platform  2 

Range  bias,  nmi  1  -2 

Azimuth  bias  0  0 

Elevation  bias  0  0 

The  results  are  illustrated  in  Fig.  2.  Approximately  ten  iterations  were  required  for  the  filter  to  accu¬ 
rately  obtain  the  induced  biases. 

These  results  verify,  under  simulation,  that  a  simple  algorithm  based  on  requiring  two  non- 
colocated  radars  to  measure  the  same  distance  between  targets  in  space  can  accurately  estimate 
range  bias  errors.  More  realistic  testing  is  required  to  verify  operation  of  the  algorithm  in  ship¬ 
board  environments.  In  the  next  section  the  results  of  testing  with  real  data  will  be  discussed. 


(18) 


(19) 


7 


PRATHER 


so 

o» 


«*•  o  i 

LO  V  "I 


Fig.  2  —  Induced  range  bias 


RESULTS  WITH  REAL  DATA 

Further  tests  were  made  on  the  range  bias  algorithm  using  detection  data  of  aircraft  targets 
of  opportunity  recorded  during  the  month  of  September  1979  by  the  SPS-39  radars  located  at  CBD 
and  at  JHU/APL.  To  test  the  algorithm,  a  specific  track  was  selected  which  was  no  less  than  16  scans 
in  duration.  This  selected  track’s  position  along  with  a  second  randomly  selected  track’s  position 
produced  the  difference  vectors  which  were  used  for  input  to  the  least-square  filter.  At  the  end  of 
each  scan,  NRL  data  were  transformed  from  the  NRL  coordinate  system  to  that  of  APL’s  for  com¬ 
parison  with  and  without  biases  added.  Bias  removal  calculations  [2]  were  performed,  and  the  fol¬ 
lowing  bias  errors  were  used  in  conjunction  with  the  range  bias  errors  estimated  by  this  algorithm: 

on  the  APL  radar  2.033°  in  azimuth 

0.119°  in  elevation 

on  the  NRL  radar  -0.691°  in  azimuth 

0.505°  in  elevation 

Applying  these  bias  errors  to  several  tracks  produced  results  that  were  very  encouraging.  Refer¬ 
ring  to  Fig.  3,  we  have  a  track  as  seen  by  both  sites  without  any  biases  removed.  Figure  4  shows  the 
results  of  removing  the  above  azimuth  and  elevation  biases  plus  the  following  range  biases: 

on  the  APL  radar  1.895  km  (1.162  mi) 

1.714  km  (1.071  mi) 


on  the  NRL  radar 
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These  range  biases  were  obtained  by  taking  the  average  on  biases  obtained  on  several  runs.  In  Figs.  3 
and  4,  as  well  as  the  others,  APL  is  located  at  the  origin  with  NRL  located  approximately  56.5  km 
(35  mil  southeast.  Figures  5  and  6,  again,  represent  a  track  as  seen  by  both  sites  without  and  with 
biases  respectively.  Figures  7  and  8  are  similar  to  the  above  but  produced  results  which  were  not  as 
accurate. 


SUMMARY 

An  algorithm  was  developed  to  measure  the  range  bias  errors  that  can  occur  in  radars.  The 
method  was  based  on  the  condition  that  two  non-colocated  radars  must  measure  the  same  distance 
between  targets.  Because  only  distances  between  points  were  being  measured  and  compared,  it  was 
not  necessary  to  know  the  radars’  location  or  orientation.  This  fact  is  advantageous,  because  the 
entire  “Gridlock”  problem  can  then  be  solved  in  independent  stages. 

In  this  study  the  elevation  bias  error  was  assumed  to  be  zero  or  previously  removed.  However, 
in  a  more  expanded  study  this  requirement  is  not  necessary  and  could  be  included  in  the  same 
estimation  process  presented.  The  linearization  of  the  equations  seemed  to  work  quite  well  and  did 
not  present  a  problem  with  the  types  of  data  studied. 

Excellent  results  were  obtained  in  removing  bias  errors  under  controlled  simulation  conditions. 
However,  when  the  recorded  data  from  the  APL  and  NRL  sites  were  used,  the  results  did  not  look 
quite  as  good.  Time  did  not  permit  a  full  assessment  of  the  discrepancies  found  in  the  real  data. 
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Appendix 

MATRIX  COEFFICIENTS 


The  matrix 


CM 


bD  , 
bb. 


is  given  by 


bDl  by/T  1  bF 

bb  '  bb  ~  2y/F  bb  ’ 

ri  ri  V  rl 


(Al) 


where 


F~rii  +ri2  -  2riiri2[sineen  sin  eC12  +  cos  eCn  cos  e£12  cos  (aZll  -  az12)j  .  (A2) 

Recalling  from  Eq.  (3)  the  models  of  ^  2’  1  >  etc  >  a°d  ^en  differentiating  Eq.  ( A2)  with 

respect  to  br  yields, 


bD 

bb 


r~^F^ri1  +2rX2“  2sineeu  sin  cei2(rll  +  r12> 

-  2coseen  cose£12  cos  (ozn  -  oz12)(r11  +r12)J  . 


(A3) 


Collecting  terms  yields 


bD 


db  ~ ^(ru  +ri2>(1' sine8u  sineCi2  "  coseCu  coseCi2  cos  (ozn  -  «12))J.  (A4) 
By  replacing  the  subscripts  in  Eq.  (A3)  we  arrive  at  [C2i] 

^2  1  r 

"7/jr  [(r21  +r22K1  "  sin  eC2l  sin  eC22  ■  cos  e221  COS  e^22  cos  (az21  ~  «*220|  (A5) 

r2  v  2  j 


bD 


bb 


'  2 
where 


F2  =r2i  +r22  "  2r2ir22[sin  e82i  sin  e?2  2  +cose£21  coseS22  cos(oz21  -  az22)J  . 
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The  matrix 


Cci  2] 


3D, 

bDx 

3D, 

3D, 

bl) , 

3D 

3  n 

dn 

3  n 

bn 

3n  ,, 

2  ‘f, 

3/i  , 

r\  1 

r\ 2 

az\  1 

1 

1 

is  a  matrix  made  up  of  position  elements  differentiated  with  respect  to  the  measurement  noises. 
Recalling  again  from  Eq.  (3)  the  models  of  r, , ,  r,  2.  1 ,  etc.,  ami  expanding  the  first  element 

yields. 

3D,  by/F  j 
b n  bn  2  JT  3 n 

rtt  rn  rn 

with  F  from  Eq.  ( A2). 

Differentiating  Eq.  (A2)  with  respect  to  n  yields 

'11 

3  D 


'  1 1 


1  1  r 

"  =  2s/F  t2 

2r,  2  cos  efl,  ,  cos  et!,  2  cos  (az, ,  -  az,  2)J  . 


3 „  2JT  l2r'  1  2r!2sin<*Vll  s'n^,2 

rt  1  ^ 


Collecting  terms  yields 
3D 


1  _.  A.  [V,  ,  -  r,  2  (sin  1  +  sin  ^'•’i 2  +  cos e^i  1  cos 2 cos  < az\  1  az\ 2OI 

3n  y/F  L  J 


Similarly  differentiating  the  remaining  elements  of  C,  2  yields 
3D, 


dn 

ri  2 

y/T  L 

3D, 

bn 

az  1 1 

s/T  [ 

3D, 

1 

3n 

az,2 

"“v* 

3D, 

-  —  fr 

bn  f 
eKn 

~  s/F  L 

1 

v^r[rllrl2  C08^ii  cos  e«,  2  sin  (az, ,  -  az,2)]  , 


1  1  r  I 

—  =  “  jp  |_rnri2  l'os  *®i  1  cos  p?i  2  sln  <a*i )  0*1 2*J  • 


=  [r,  ]  rx  2  [sin  et!, ,  cos  eC,  2  cos  (az,  ,  az,  2)  -  cos  ct,  ,  sin  ef,  2jJ  . 


(Ad) 


( A7) 


“  =  W[ri  2  “  ri  1  (sin  <’ci  1  sin  2  +  coseV, ,  coseC,  2  cos  (az,,  -  az,2))J  ,  (A8) 


(A9) 


(A10) 


(All) 
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I—  rrnr12[coseBn  sin  e«12cos  («,,  -  ^12)-  sin  eCn  coseC,2]J 
dn  q  n/P  L 


( A12) 


By  replacing  the  subscripts  in  Eqs.  (A6)  through  (All),  we  arrive  at  the  elements  of  C22>  with  F2 
defined  by  Eq.  ( A5): 

&D  "1 

- 2_=  J— [r  -  r22[sineC2,  sin  e£22  +cosee21  cos  e822  cos  (02 21  -  «22)]J  ’  (A13) 

3nr21  L 

— — =  —y=  ["r22  -  r21  [sineC21  sin  e222  +  cos  eC21  cosetM  cos  (az.n  -  «22)]]  -  (A14> 

3nr00  VF2  L 


3Do  1  -1 

- =  -^[r21r22  cose821  COS  eC22  Sin  “21  J  > 

'“*21  ^  2 

y 

— —  =  -  -7==-  fr2ir22  coseC21  cos  efi22  sin  <w22]  . 

-  -^rrr2ir22  [sin  e«21  coseC22  cos  ( az2l  -  a*2 2)  "  coscC21  sm  ^22]]  ’ 

M.,  v*2  L 


(A15) 


( A16) 


(A17) 


M>2  1 


=  -ip £  r21r22  [cos  e821  sin  e622  cos  <«2l  -  <*22)  -  sin  e$>21  cos  el>22]j  .  (A18) 
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